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Abstract: The uncertainty of nuclide libraries in the analysis of the gamma spectra of 
low- and intermediate-level radioactive waste (LILW) using existing methods 
produces unstable results. To address this problem, a novel spectral analysis method is 
proposed in this study. In this method, overlapping peaks are located using a 
continuous wavelet transform. An improved quadratic convolution method is 
proposed to calculate the widths of the peaks and establish a fourth-order filter model 
to estimate the Compton edge baseline with the overlapping peaks. Combined with 
the adaptive sensitive nonlinear iterative peak, this method can effectively subtracts 
the background. Finally, a function describing the peak shape as a filter is used to 
deconvolve the energy spectrum to achieve accurate qualitative and quantitative 
analyses of the nuclide without the aid of a nuclide library. Gamma spectrum 
acquisition experiments for standard point sources of Cs-137 and Eu-152, a 
segmented gamma scanning experiment for a 200 L standard drum, and a Monte 
Carlo simulation experiment for triple overlapping peaks using the closest energy of 
three typical LILW nuclides (Sb-125, Sb-124, and Cs-134) are conducted. The results 
of the experiments indicate that (1) the novel method and gamma vision (GV) with an 
accurate nuclide library have the same spectral analysis capability, and the peak area 
calculation error is less than 4%; (2) compared with the GV, the analysis results of the 
novel method are more stable; (3) the novel method can be applied to the activity 
measurement of LILW, and the error of the activity reconstruction at the equivalent 
radius is 2.4%; and (4) The proposed novel method can quantitatively analyze all 
nuclides in LILW without a nuclide library. This novel method can improve the 
accuracy and precision of LILW measurements, provide key technical support for the 
reasonable disposal of LILW, and ensure the safety of humans and the environment. 
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1. Introduction 


A 1000 MWe nuclear power unit will produce approximately 13,000 m° of low- 
level and intermediate-level radioactive waste (LILW) !! after 60 years of operation. 
After solidification or compression, this waste is prepared in 200 L and 400 L drums 
7], with a total volume of approximately 3000 m°. An International Atomic Energy 
Agency (IAEA) technical report P! stated that to ensure safety within the entire LILW 
management life cycle, the LILW must be characterized. At the same time, according 
to the IAEA safety standard l], the disposal depths of LILW for different activity 
ranges differ, and thus quantitative analysis of LILW is necessary. Because destructive 
methods and time-consuming drilling are necessary to reopen a container M“, a 
nondestructive assay (NDA) is generally the best choice. Some important LILW 
nuclides include alpha-emitting nuclides (most actinides and transuranic nuclides) and 
beta-emitting nuclides (C-14, Ni-63), which cannot be used for NDA. Thus, these 
nuclides are called difficult-to-measure (DTM). According to the International 
Organization for Standardization l, by performing gamma scanning on easy-to- 
measure nuclides (also called key nuclides) such as Cs-137 and Co-60, it is possible 
to establish a scaling factor between key nuclides and DTM and calculate the activity 
of the DTM. Therefore, gamma-ray scanning is required to characterize LILW. 

According to the American Society for Testing and Materials (ASTM) Standard 
Guide '*!, NDA methods include segmented gamma scanning (SGS) '7! and 
tomographic gamma scanning (TGS). TGS is particularly well-suited for items with 
heterogeneous matrix materials and non-uniform radioisotope distributions. Some 
studies have been conducted on SGS and TGS [88], For both SGS and TGS, it is 
necessary to use gamma ray spectroscopy to conduct qualitative and quantitative 
analyses of the LILW gamma spectrum. Owing to the diverse range of nuclides in 


[19,20] 


LILW and the variations between waste streams , Identifying nuclides can be 


challenging. Therefore, high-resolution high purity germanium (HPGe) gamma ray 


spectrometry is the optimal method for LILW measurement 1-251, 
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The spectrum measured by an HPGe detector must be analyzed to obtain the 
counts of gamma rays with different energies. The main steps of the spectrum analysis 
method °! include smoothing, background subtraction 7, peak searching, and 
calculation of the peak area. Nuclide identification relies on the process of peak 
searching. Currently, the most mature peak-searching method is the derivative method 
proposed by Mariscotti °*!. GammaVision (GV) ! and Genie 2000 °°! have adopted 
this method. Because of the limited peak-searching ability of the derivative method, it 
is necessary to establish a nuclide library to assist in identifying overlapping peaks 
composed of two full-energy peaks that are very close to each other. This method 
assumes that all gamma rays in the nuclide library are present in the spectrum, and 
thus it attempts to fit the peaks for every gamma energy listed in the library. Therefore, 
for the same peak, different results are often obtained with different nuclide libraries. 
For example, common LILW nuclides include Cs-137 emitting 662 keV and Ag-110m 
emitting 658 keV. The nuclide composition of the LILW before measurement is not 
accurately known. For the same full-energy peak of 662 keV, if only 662 keV is 
defined in the nuclide library, the quantitative analysis result is accurate. However, if 
both 658 keV and 662 keV are provided in the nuclide library, the report will provide 
a full-energy peak count of 658 keV. This results in misidentification of the nuclide. 
At the same time, the calculated 662 keV peak area decreases, resulting in inaccurate 
activity calculations, which is a potential safety hazard for the disposal of LILW. The 
IAEA report 7?! also noted that this method requires the availability of good nuclear 
data libraries. Other peak-finding methods include the IF function method PH, zero- 
area fold-in method P”, Gaussian product function method 6°31, covariance method 4, 
continuous wavelet transform (CWT) method 05I, and neural network method 2°. 
Compared with the derivative method, these methods are an improvement; however, 
they still require auxiliary identification of the nuclide library. 

In this study, to solve the problem of uncertainty in the nuclide library causing 
instability in the analysis results for LILW spectra, a novel spectrum analysis method 


is proposed to reduce this uncertainty and improve accuracy. For each step in the 
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spectral analysis process, the overlapping peak decomposition process does not rely 
on the assistance of the nuclide library. This study improves the accuracy of LILW 


measurements and provides key technical support for the safe disposal of LILW. 


2. Materials and Methods 


2.1. Analysis method 
2.1.1. Traditional spectral analysis method 

In traditional spectral analysis methods, there are problems with each of the basic 
steps. When smoothing the spectrum according to a five-point smoothing algorithm 
[37], the processing causes the maximum value of the peak to decrease each time, and 
the counts on both sides of the bottom of the peak increase. Thus, the calculation of 
the peak area is inaccurate. At the same time, it is difficult to identify the nuclide 
because of the changes in the peak shape. 

According to the GammaVision user manual °7!, background subtraction methods 
include the straight-line, stepped, and parabolic background methods "*!. These 
methods provide only rough approximations of the background, thereby increasing 
the peak area calculation error. At the same time, Miroslav °! showed that the 
parabolic background had very poor robustness B®, 

The peak search method was proposed by Mariscotti ?°!. This method has a 
limited ability to separate overlapping peaks '7!, resulting in inaccurate nuclide 
identification; furthermore, it is highly susceptible to noise. Thus, it is necessary to 
smooth the spectrum multiple times before searching for the peak, which leads to 
changes in the peak shape and increases the error in the peak area calculation. 

The final step is the quantitative analysis of overlapping peaks. In traditional 
methods, a preliminary library-based peak search of the spectrum is performed. The 
analysis assumes that all of the gamma rays listed in the library are present in the 
spectrum; therefore, an attempt is made to fit a peak at every gamma energy listed in 


the library. Therefore, as reported in Ref. '*°!, if the nuclide library is inaccurate, it will 
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increase the error in the spectral analysis results. 


2.1.2. Spectral analysis method proposed in this study 

According to Knoll | and the IEC international standard '7!, for a single 
overlapping peak composed of multiple full-energy peaks with similar energies, the 
only method to conduct an analysis is through deconvolution. The deconvolution 
method can be used to establish the relationship between the incident gamma ray 
spectrum and the energy deposition spectrum using the response function of the 


detector, as follows !*8!: 

y(t) = i x(t)h(t — t)dt + n(t), (1) 
where x(t) is the incident spectrum, y(t) is the energy deposition spectrum, h(t) is 
the response function, t is the channel, and n(t) is the system noise. Therefore, the 
core concept of this method is to use the peak shape function, h(t), as a filter, and 
then to obtain the incident spectrum through deconvolution to realize the 
decomposition of overlapping peaks. 

For h(t), because of the different attenuations of the medium during gamma ray 
measurement under different conditions, it is difficult to establish the full-spectrum 
response under such complex and variable conditions through experiments or Monte 
Carlo simulations. Because the full-energy peak areas of the different nuclides are 
regarded as regions of interest in this study, it was only necessary to determine the 
full-energy peak response. The full-energy peak of the HPGe detector can be 
represented by a Gaussian function. Theoretically, according to the full width at half 
maximum (FWHM) of the full-energy peak, the response can then be obtained via 
calculations. Therefore, before the spectral analysis, an experimental method was used 
to calibrate the FWHM. The FWHM is different for different energies; therefore, the 
following equation is used for parameter fitting 14: 

FWHM = a+ byE +c x E2, @) 
where E is the gamma energy. The HPGe detector used in this study was calibrated 
using the experimental method. When the energy unit is MeV, the values of a, b, and c 


are 0.00068, 0.00062, and 1.07, respectively. 
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For a discrete system, h(t) can be expressed as a matrix l5]; 


hA) 0 0; 0 
hy(2)h2,(1)_ 0 : 
hı (3)h2(2)h3(1): 

h2(3)h3(2): 
hy(A) 7 iG): hhl |, n 
O M(B) i © Mm(2) 
a N =) 


O 0 O È hn) 
where each column of the matrix is the probability density function of the full-energy 
peak response, and the subscripts of h indicate the gamma rays emitted from that 
channel. Because the FWHMs of the different full-energy peaks are different, the 
values of a—n represent the peak widths of different full-energy peaks. 

The solution to Eq. (1) is an ill-posed problem. According to Miroslav "5, 
iterative regularization has strong positive constraints and noise-suppression 
capabilities. Currently, the most effective method for gamma spectrum analysis is the 
boosted traditional iterative algorithm 7). In this study, a boosted maximum- 
likelihood expectation-maximization algorithm '4*! was used. 

Through this deconvolution process, the deposition spectrum of gamma rays with 
different energies can be converted into the incident 6 function spectrum; thus, 
nuclide identification can be easily performed. Its height corresponds to the area of 
the full-energy peak. However, in the realization of this process, another problem 
must be solved, i.e., the background must be effectively subtracted from the original 
spectrum to obtain a “net spectrum” containing only the full-energy peak. 

Background subtraction methods include the linear method °!, approximate 
physical models for the background, the polynomial method "°l, digital filtering 1%, 
peak clipping '*!, fast Fourier transform 5", and the sensitive nonlinear iterative peak 
(SNIP) method?! Among these methods, the most effective background subtraction 
method is the improved SNIP method, i.e., adaptive SNIP B°, 

The basic concept of the adaptive SNIP method is to divide the energy spectrum 
into two regions: the region where the peak is located, and the region where there is 


no peak. For the peak area, the SNIP filter is used for background subtraction, and for 
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the non-peak area, all of the gamma counts are regarded as the background. Therefore, 
the problem is to determine how to calculate the peak position, peak width, and 
selection of the clipping filter. 

This study proposes a calculation method for peak position identification based 
on CWT. The wavelet transform of a signal x is given as follows $5: 

T(E,s) =f. x(t)We,(tdt, @) 
where E is the gamma ray energy, s is the scale of the wavelet, Yg s(t) is the mother 
wavelet, and T(E,s) is the scalogram. The scale s represents the stretching and 
shrinking of the wavelet, as expressed in the following equation 131: 

ve s(t) = =o (S), ©) 
The correct choice of wavelet is the key to obtaining the results of the spectrum 
analysis. As the full-energy peak of the spectrum is a Gaussian function, a Mexican 
hat wavelet (Mexh) is selected because its shape is similar to that of a Gaussian peak. 


The definition of Mexh is given as follows 4): 
t2 


Puexht) = uae — ©) e 20%, (6) 


301 


By analyzing the scalogram, the position information of the peak can be obtained. 
The next step is to calculate the widths of the peaks based on their positions. 

The peak width can be obtained directly from the scalogram obtained with Eq. (4), 
but because the peak widths of different scales are different, this method has large 
errors [55], The peak width can be calculated using a Gaussian product function °?!, but 
it is highly susceptible to noise interference. The width can also be determined from 


[56], This approach is also difficult to use 


the shapes of the peaks using function fitting 
because the number of singlets comprising the overlapping peaks is unknown. In Refs. 
[55] and 556l, a method for calculating the peak width based on the second convolution 
of the Gaussian first derivative was proposed. 

The expression for the full-energy peak is as follows ©7!: 


2 
= 7 
G(x) = Aexp (- 9 ), y) 


where x is the channel, A is the count of the full-energy peak maxima, a is the peak 


7/31 


position, and ø is the standard deviation of the full-energy peak. 
The calculation result for the second convolution of the Gaussian first-order 


derivative of Eq. (7) is expressed as follows 4): 


__ 2n6*o((x—a)?-0?—26") _ (x—a)* (8) 
CC) = aap OP (anan) 


where ô represents the width of the filter. 


Using Eq. (8), the crossing point of the zero value can be determined as follows 


[54]. 


o = J (Xo — a)? — 28, ©) 


where a a is the peak position calculated with the CWT, and xq is the intersection 
point with the x-axis. Then, [a — 30,a+ 3a] is the width of the area in which the 
peak is located. 

However, other studies have shown that this method is not suitable for calculating 


39]. Therefore, in 


overlapping peak widths that contain too many full-energy peaks ! 
this study, we propose an improved method. 

For Eq. (8), the calculation of the width of the peak requires the determination of 
the parameter o, while the parameter 6 is also unknown. Thus, a method is required to 


determine parameters o and 6. For C(x) in Eq. (8), when ae = 0, it is easy to obtain 


ô = o, and this is the maximum value of C(x). Therefore, the absolute value of C(x) 
is first computed, and then a matrix M(x, 6) is built. Its horizontal axis, x, represents 
different channels. The vertical axis of M(x, ô) indicates that for each x, the value of 
C(x) is calculated for different 5. In other words, for each x, the calculation is 
performed by traversing the possible values of 6 (for example, 0.1-100) and then 
calculating a C(x) curve. This is carried out because at the maximum value of C(x), 
ô = 0, the value of o can be obtained based on the value of 6, and the peak width can 
then be effectively calculated. In the application, the horizontal axis, x, only considers 
the peak position information calculated based on the CWT. 

Next, the background of the peak area must be removed. For the SNIP method, 


first, to compress the dynamic range of each channel count in the spectrum to speed 
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up the iteration, a linear least-squares operator is applied to the spectrum data !*7!: 

v(i) = log flog| ya+1it 1| + 1}, (10) 
where i is the channel, y(i) is the count of each channel, and v(i) is the working 
vector. 

For the count v(i) of each channel, the parameter m is selected, and then w = 2m 
+ 1. The following equation is then used to iteratively calculate v, (i), v2 (i), ..., Vm 
[52]. 
vp (i): = min {vp 0,5 [vpa(E + p) + v1 - PII}, 
where p E [1, m], i € [p, size,], and size, is the size of spectrum y. For parameter w, 


when the value of w is too small, the background estimate will be too large, and the 
net peak area will be too small. If the value of w is too large, the background estimate 
will be small, and the net peak area will be large. Theoretically, the background 
should be removed as much as possible while preserving the net peak areas. The real 
width of an object (peak, doublet, or multiplet) that should be preserved is the best 
value of w 53], In the proposed method, the width of each peak is calculated based on 
the second convolution of the Gaussian first derivative, and the clipping width, w, of 
each peak is obtained; that is, the value of w is dynamic. Then, the SNIP algorithm is 
utilized. Finally, a linear least-squares inverse operator operation is performed on 
Vm(i) to obtain the background spectrum. 

For the adaptive SNIP D°], the spectrum is divided into two areas: the non-peak 
and peak areas. In the non-peak region, all of the gamma counts are the background. 
In the peak region, the background is subtracted using the traditional SNIP method. 
Eq. (11) represents a second-order filter that can effectively remove the linear 
background. To solve the problem of the inability to subtract the background at the 
Compton edge for the full-energy peak as described in Section 2.1.1, a novel method 
is proposed. 

Because the Compton edge can be approximated by a cubic function to some 


extent, and a cubic function can be removed by a fourth-order clipping filter, a fourth- 
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order clipping filter is proposed for background subtraction in the peak area. The 
parameters of the equation were derived from Pascal’s triangle ©?! formula. The 


fourth-order clipping filter is defined as follows: 


Vi = Vp- Ù 
V, = Vp—1(i—p)+vp-1(i+p) 
i 2 (12) 
—Vpy—1 ip) +4 Vp -1 (iE) H4vp-1 (i+4) —Vp_1(i+p)’ 
6 
v(i) = min [V1, max (V2,V3)] 


V3= 


where the definitions of these parameters are identical to those in Eq. (11). 

Through the peak searching and peak width calculations described above, the 
spectrum can be divided into peak and non-peak areas, which are combined with a 
fourth-order clipping filter in the peak area. Thus, effective subtraction of the 
background is realized and the net spectrum is obtained. 

During the smoothing process in the traditional method, the peak count will 
decrease. This study uses a low-pass filter 7! to smooth the spectrum and effectively 
solve this problem. 

2.2. Experimental condition 
2.2.1. Experimental configuration 

The HPGe detector was an Integrated Cryocooling System produced by ORTEC. 
Its relative efficiency at 1.33 MeV is 30%. The detector crystal has a diameter of 51.8 
mm and length of 60.9 mm. The cold finger has a diameter of 8.7 mm and a length of 
38.2 mm. The dead layer at the front end of the crystal has a thickness of 0.03 mm. 
The DSPEC-50 multi-channel analyzer was manufactured by ORTEC. 

GV 6.08 spectral analysis software was used. The analysis parameters of GV 
included the calibration (.CLB) file, peak background correction (.PBC) file, and 
nuclide library (.LIB) file; the analysis program was WAN32. The standard area of 
the peak was calculated by selecting the region of interest (ROI) in the GV; the peak 
width was three FWHM. 


The experiments were conducted using the Tomographic Gamma System (SJTU- 
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TGS), which was independently developed by Shanghai Jiao Tong University, as 
shown in Fig. l(a). The detector was equipped with a lead collimator (liftable). The 
exterior dimensions of the collimator and collimation window were 40 cm (L) x 36 
cm (W)x 25 cm (H)x and 25 cm (L) x 20 cm (W)x 9 cm (H), respectively, and the 
thickness of the lead was 8 cm. The distance between the detector surface and drum 
center was 74 cm. The system had a liftable transmission source storage device. The 
waste drum could be moved and rotated. Therefore, the measurement system had 
detection functions of SGS and TGS and could also meet the requirements of 200 L 


and 400 L waste drums at the same time. 


(a) (b) 


Fig. 1 (a) SJTU-TGS System; (b) Canberra SGS System 


2.2.2. Source test conditions 

To verify the accuracy of the spectral analysis of the different nuclides, the 
radioactive sources used were standard point sources: monoenergetic gamma ray 
nuclide Cs-137 !! (7.86E+05 Bq) and multi-energy gamma-ray nuclide Eu-152 [61 
(2.74E+06 Bq). 

To obtain the spectrum of the radioactive source attenuated by different densities, 
the radioactive source was placed in a 200 L drum filled with different materials for 
measurement. The materials used in the drum were air (the empty drum), water (1 
g/cm), and sand (1.6 g/cm’). 

To explore the influence of different nuclide libraries on the GV analysis results, 


multiple verification experiments were conducted on 200 L drums. The radioactive 
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source, Cs-137, was placed at different positions in the drum, as shown in Fig. 2, and 


the drum was not filled with material. The time for each measurement was 10 min. 


HPGe detector 


pon | 200 L drum 


Fig. 2 Schematic diagram of the locations of the radioactive source 


2.2.3. SGS experiment 

The purpose of spectral analysis is to measure the activity of LILW. Therefore, a 
200 L standard drum was employed. The material in the drum was water. Cs-137 was 
placed at radii of 0, 6.25, 11.25, 16.25, 21.25, and 26.25 cm, and the height was 42.5 
cm (half the height of the drum). SGS 1-161] was used for activity reconstruction based 
on the SJTU-TGS system. 

The same standard drum was used for activity reconstruction with the Canberra 
SGS System [5], as shown in Fig. 1(b). The aim of this study was to compare and 
analyze the activity measurement results of the two systems. The spectroscopy 
software in the Canberra SGS System is Genie 2000 B%, and the software for 
measurement activity is NDA2000 '!. To control for the variables, the SJTU-TGS 
System adopts the SGS method described in NDA2000 (refer to Refs. I! and '°!!). The 
efficiency calibration method of the SJTU-TGS system adopts the method described 
by Qian |'7!, which meets the ASTM standards [%1]. Because the material was water, 
this was an easy task for efficient calibration. It was assumed that any calibration 
method would provide relatively consistent results. During the gamma scanning, the 
drum was vertically divided into eight layers, and the lifetime of each layer was 1 min. 
Considering the dead time and moving time of the detector, the measurement time for 
a drum was approximately 10 min, which was within the measurement time range for 
200 L drums in nuclear power plants. The only difference is that the SJTU-TGS 
System adopted the novel spectral analysis method proposed in this study. For the 
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Canberra SGS System according to Genie 2000, its method was very close to GV, and 


thus it represented the traditional spectrum analysis method. 


3. Experiments and Discussion 


The verification experiments for the spectral analysis method included the 


following parts: a preliminary verification experiment for a single nuclide, instability 


experiments for the nuclide library, activity measurement comparison experiments for 


the standard drum, and overlapping peak simulation spectrum experiments. 


3.1. Preliminary verification experiment for a single nuclide 


3.1.1. 


Spectra of monoenergetic gamma ray nuclides 
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Fig. 3 (a) Original spectra of Cs- 
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MM observed spectrum 
1.0 mm after decomposition 
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(c) 
Fig. 4 Analysis results of the Cs-137 spectrum with different methods: (a)-(c) spectra attenuated by different 


materials (air, water, and sand). Left: Standard peak area of ROI; middle: GV; right: Novel method 


Table 1 Analysis results for the Cs-137 spectrum attenuated by different materials with different methods 


Standard 
Material Method peak area Calculated peak Error 
area 
of ROI 
GV 55676 2.4% 
Air 54226 
Novel method 52895 —2.71% 
GV 1091 0.65% 
Water 1084 
Novel method 1065 —1.75% 
GV 1004 2.66% 
Sand 978 
Novel method 982 0.41% 


The original spectra of the Cs-137 nuclide without attenuation and those 
attenuated by water and sand are shown in Fig. 3(a); the normalized spectra are shown 


in Fig. 3 (b). The Compton plateau in the attenuated spectrum is significantly higher 
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than that in the unattenuated spectrum. Three spectra were analyzed using the 
proposed method. For comparison, GV was used to establish an accurate nuclide 
library and the nuclide analysis results were obtained. The spectral analysis results 
obtained using different methods are shown in Fig. 4, and the quantitative results are 
presented in Table 1. 

As shown in Fig. 4, through the deconvolution process, the proposed method 
converted the number of full-energy peak channels into one channel, making it very 
easy to identify nuclides. The results in Table 1 indicate that because of the simple 
task of single-peak nuclide identification, the results were similar to those of GV, and 


the errors for both methods were less than 3%. 


3.1.2. Spectra of multi-energy gamma ray nuclides 
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Fig. 5 (a) Original spectrum of Eu-152; (b) Spectrum analysis results for Eu-152 with the novel method. 
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The percentages represent gamma ray emission probabilities. 


Table 2 Comparison of the spectral analysis results for Eu-152 with different methods 
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Energy Standard | Calculated peak | Error (GV) Calculated peak Error (novel 
peak area area (GV) area (novel method) method) 
121.78 47650 48026 0.79% 49262 3.38% 
244.70 10908 10844 -0.59% 11115 1.90% 
344.28 30823 30083 -2.40% 31042 0.71% 
411.11 2283 2703 18.40% 2401 5.17% 
443.97 2873 2987 3.97% 3111 8.28% 
778.90 7794 8057 3.37% 8047 3.25% 
867.37 2177 2456 12.82% 2214 1.70% 
964.08 7377 7497 1.63% 7811 5.88% 
1085.87 4760 5235 9.98% 4953 4.05% 
1112.07 6369 6704 5.26% 6589 3.45% 
1408.00 8122 8163 0.50% 8528 5.00% 
Mean 5.43% 3.89% 
error 


Because the nuclides in LILW generally emit multi-energy gamma rays, the Eu- 
152 spectrum was used for analysis and verification. The original spectrum of Eu-152 
is shown in Fig. 5 (a). Gamma rays for eleven energies with the largest peak areas 
were selected for analysis. Because the energy resolution of HPGe is very good, the 
overall spectral lines appeared very clear. However, in the zoomed-in details, such as 
the 778.90 keV full-energy peak shown in Fig. 5 (a), the peak width was 
approximately 20 channels. The results of the spectral analysis using the novel 
method are shown in Fig. 5 (b). It can be observed that the gamma rays occupied only 
one channel, as seen in the zoomed-in details of the 778.90 keV full-energy peak. 
Thus, the accurate identification of nuclides could be realized. 

The results of the quantitative analysis of the novel method and GV are compared 
in Table 2. For the calculation errors of the full-energy peaks with different energies, 
the proposed method was more accurate, with a maximum error of 8.28%, whereas 
the maximum error of GV was 18.40%. This is because in the energy spectrum 
analysis process of GV, spectral smoothing is first performed, leading to some 
changes in the peak shape, and the peak area is calculated by fitting these peak shapes. 
Therefore for peaks with large statistical fluctuations, particularly those with 
relatively small peak areas, significant errors can be introduced during the smoothing 
and fitting processes. However, in the novel method, the full-energy peak is directly 
used to calculate the peak area during the process of deconvolution; thus, there is 
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almost no significant deviation in the calculation of the peak area. From the 
perspective of the mean error (average error from all the peaks), GV has an error of 


5.43%, whereas the novel method has an error of 3.89%. 


3.2. Nuclide library instability experiments 

These GV analysis results were based on an accurate nuclide library. However, in 
practical applications, it is generally difficult to make an accurate judgment regarding 
the suitability of the nuclide library. Therefore, for the same Cs-137 (662 keV) 
spectrum attenuated by sand, different nuclide libraries were used to analyze the GV 
results. Analytical results are also provided for this novel method, which does not 
require the establishment of a nuclide library. Because Ag-110m is the most common 
corrosion activation product in LILW and its gamma emission of 658 keV is close to 
that of Cs-137, the nuclide library settings were as follows: Cs-137, Cs-137+Ag-110m, 


and Ag-110m. The GV analysis results are presented in Fig. 6 and Table 3. 
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(a) (b) (c) (d) 
Fig. 6 (a) Cs-137 original spectrum; (b) GV results based on the nuclide library for Cs-137; (c) GV results 


based on the nuclide library for Cs-137+Ag-110m; (d) GV results based on the nuclide library for Ag-110m. 


Table 3 Comparison of spectral analysis results with different nuclide libraries 


Calculated 

Number Nuclide library | Standard Calculated Error peak area 
(GV) peak area | peak area (GV) | (GV) (novel 

eS 


pots 137 Ag- 
S- 8 32.3% 0.4% 
110m 
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As shown in Fig. 6, for the Cs-137 spectrum, if only Cs-137 was included in the 
nuclide library, the analysis result was accurate and the peak area error was only 2.7%. 
However, when Ag-110m was added to the nuclide library, a peak area of Ag-110m 
was obtained based on the spectral analysis mechanism; in other words, Gaussian 
functions with energies of 658 and 662 keV fit the full-energy peak of Cs-137. 
Therefore, the areas of the full-energy peaks of the two energies could be obtained, 
and thus the error of the peak area of Cs-137 increased to 32.3%. In the absence of 
Cs-137 in the nuclide library, the result did not provide a count of Cs-137; in other 
words, different nuclide libraries yielded completely different analytical results. By 
assuming the presence of a peak a priori, GV was “forced” to calculate the properties 
of the peak. This finding was consistent with the conclusions of Santoro '°!, Because 
the novel method does not require a nuclide library, the error in the analysis result was 
stable at 0.4%. 

To further explore the influence of different nuclide libraries on the GV analysis 
results, verification experiments were carried out on 200 L drums, as shown in Fig. 2. 
For the nuclide library setting of the GV, in addition to Cs-137 (662 keV) and Cs- 
137+Ag-110m (662 keV + 628 keV), an artificial nuclide library was used with a 
composition of 662, 661, 660, 659, and 658 keV. The results of the spectral analyses 
with different methods are shown in Fig. 7, where the x-axis represents the 
identification number of the radiation source at different positions inside the 200 L 
drum. The analysis results of the novel method and the GV based on an accurate 
nuclide library were similar. The average relative error of the novel method was 3.3%, 
and that of GV was 3.4%. However, when the GV nuclide library consisted of two 
gamma energies, the 662 keV peak was fitted into two peaks, and thus the average 
error of the 662 keV peak was 38.3%. When the nuclide library consisted of five 
gamma energies, the 662 keV peak had the worst average error of 48.3% because 662 


keV was “split” into five overlapping peaks. 
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Fig. 7 Comparison of the results of different methods 


3.3. Activity reconstruction for a 200 L drum 
For the 200 L standard drum, the SGS activity reconstruction results using the 


SJTU-TGS and Canberra SGS systems are shown in Fig. 8. 
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Fig. 8 Comparison of the activity reconstruction errors for 200 L drums with different systems 


For the two systems, except for the spectral analysis method, the other methods 
are basically the same. Therefore, theoretically, the activity reconstruction error 
primarily reflects the error of the spectral analysis method. As shown in Fig. 8, when 
the radius was relatively small, the measured value was low. As the radius increased, 
the measured value also increased. Because the SGS method was aimed at LILW in 


which the matrix and nuclides are uniformly distributed, the object of efficiency 
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calibration was a uniform volume source. The Cs-137 source at different positions in 
this experiment was categorized as a situation in which the distribution of nuclides 
was extremely uneven, which does not satisfy the basic conditions of SGS. Thus, 
relatively large measurement errors inevitably occurred. For non-uniform cases, SGS 
can reach a maximum measurement error of more than 500% [64]. Therefore, when the 
radius of the source is small, the value of the efficiency calibration is too large and the 
measured value is too small. In theory, for the SGS method, there is generally an 
equivalent radius at the position of 0.7—0.8 R (R is the radius of the drum) t5161, and 
the actual efficiency value at this position is close to the efficiency of the uniform 
source value. As a result, better measurement results can theoretically be obtained. 
Therefore, for the SJTU-TGS System, when the radius was 21.25 cm, it was close to 
the equivalent radius, and the spectral analysis results were relatively accurate. 
Therefore, the error was only 2.4%, whereas for the Canberra SGS System, the 
smallest error was 12.9%. To some extent, this reflects the accuracy of the spectrum 
analysis. 
3.4. Analysis of overlapping peaks 
3.4.1. Verification of the accuracy of the Monte Carlo simulation method 

The types of nuclides in LILW are complex, and it is difficult to obtain different 
overlapping peaks experimentally to verify the spectral analysis method. Therefore, a 
Monte Carlo simulation was used in this study to simulate the process of measuring 
the radioactive source with an HPGe detector. A geometric model was constructed 
using the Geant4 tool !6564, According to the certificate of calibration (SRS number: 
116230 and product code: 8403-GF-PNT) corresponding to the HPGe detector 
provided by ORTEC, the hybrid source (a point source in tape in a 2-inch aluminum 
ring) was placed horizontally 30 cm from the front end of the detector, and the 
detection efficiency curve was obtained experimentally. A comparison between the 
Monte Carlo simulation efficiency [1-7] and the experimental efficiency l! is shown 
in Fig. 9(a). Except for the relative error of the full-energy peak efficiency at an 
extremely low energy (46.52 keV), which was 9.89%, all of other calculation errors 
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were less than 5%. 
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Fig. 9 (a) Comparison between the Monte Carlo simulation and experimental efficiencies; (b) Comparison of 


The Cs-137 spectrum experiment was conducted according to the measurement 
standard requirements of ANSI l? and IEEE !!. The collected full-energy peak count 
exceeded 20,000. The gamma count of the environmental background within the same 
measurement time was subtracted using the channel-by-channel method to obtain the 
net spectrum. A comparison between the normalized spectrum of the experiment and 
the Monte Carlo simulation is shown in Fig. 9(b). In addition to certain statistical 


fluctuations in the experimental spectrum, the full-energy peak and Compton plateau 


Cs-137 spectra with different methods 
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and (b). 


were in good agreement, verifying the accuracy of the Monte Carlo simulations. 


3.4.2. Spectral analysis of key LILW nuclide overlapping peaks 
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Fig. 10 Analysis results of overlapping peaks with different peak area ratios: (a)-(c) represent peak area 
ratios of 5:1:10, 10:1:20, and 20:1:40, respectively. Left: original spectra; middle: spectral analysis results with the 


novel method; right: scalogram of CWT peak searching 
Table 4 Spectral analysis results with different overlapping peaks 


True value Calculated value 


Number Energy Energy Peak area 


Nuclide Peak area Peak area 
(keV) (keV) error (%) 


Sb-125 601 109278 601 108969 -0.28% 
(a) Sb-124 603 21920 603 22026 0.48% 
Cs-134 605 217095 605 216716 -0.17% 


Based on the key nuclide types of the LILW provided by the Qinshan Nuclear 
Power Plant and through the comparison and analysis of nuclides and their 
characteristic energies, three nuclides with the closest characteristic energies were 
selected: Sb-125 (601 keV), Sb-124 (603 keV), and Cs-134 (605 keV). 

Different peak area ratios were set for the compositions of the three overlapping 
peaks. Because the overlapping peaks are difficult to decompose, it is easier to 
identify the nuclide as a full-energy peak located in the middle channel of the 


overlapping peaks. Therefore, if the peak area at the middle energy is lower than that 
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on both sides, thus forming a weak peak, it is easier to elucidate the performance of 
the spectral analysis methods. Therefore, three overlapping peaks were set with peak 
area ratios of 5:1:10, 10:1:20, and 20:1:40 for 601, 603, and 605 keV, respectively. 
For overlapping peaks with larger area ratios, because the activity calculated with this 
weak peak does not contribute significantly to the total activity of LILW and does not 
affect its final disposal, it was not necessary to continue studying these cases. At the 
same time, to illustrate the difficulty of identifying overlapping peaks, the CWT with 
the strongest performance in the traditional peak searching method was used for 
comparison. The results of the spectral analysis are shown in Fig. 10, and the 
quantitative results are summarized in Table 4. 

As shown in Fig. 10(a), through the precise deconvolution process, the novel 
method decomposed the overlapping peaks into three incident spectra whose shape 
was a ô function, and the positions of the three full-energy peaks were accurately 
identified without energy deviation. The maximum area calculation error of the strong 
peak was -0.28%, and the error of the weak peak was 0.48%. CWT was used to 
analyze the overlapping peaks, and three full-energy peaks were observed in the 
scalogram. However, one weak peak was located in the middle of the overlapping 
peaks; therefore, the weak peak was more difficult to identify. As shown in Fig. 10(b), 
there was no energy deviation in the positioning of the three full-energy peaks. The 
area maximum calculation error for the strong peak was -0.28% at most, and that for 
the weak peak was 1.72%. 

As shown in Fig. 10(c), under the condition of a peak area ratio of 20:1:40, by 
zooming in on the details in the original spectrum, because the area ratio of the three 
full-energy peaks was very different, it was determined that the 603 keV full-energy 
peak in the middle was completely covered by the high-energy peaks on both sides. 
The 603 keV full-energy peak was submerged. In the CWT scalogram, the weak peak 
at 603 keV could not be identified, and only two strong peaks at 601 and 605 keV 
could be identified. With the novel method, accurate nuclide identification could still 


be achieved for overlapping peaks in such extreme cases. The largest area calculation 
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error for the strong peak was —0.3%, whereas that for the weak peak was 3.34%. 

As shown in Fig. 10, for overlapping peaks under some extreme conditions, it 
was difficult for CWT to achieve an accurate nuclide analysis. CWT in the novel 
method was used to determine the positions of the overlapping peaks and perform 
subsequent peak width IAEA calculations. At this point, it did not matter which 


nuclides constituted the overlapping peaks. 


4. Conclusion 


In LILW spectrum analysis using existing methods, uncertainty of the nuclide 
library causes instability in the analysis results and reduces the accuracy of activity 
reconstruction. To address this problem, this study proposes a novel spectral analysis 
method. Through standard source experiments, standard drum experiments, and 
Monte Carlo simulation experiments, the following conclusions are drawn: 

1. The novel method and GV with an accurate nuclide library exhibited similar 
spectrum analysis capabilities, and the analysis errors were less than 3% for 
Cs-137 and 4% for Eu-152. 

2. The novel method does not require a nuclide library, and the analysis results 
were more stable than those of GV, which showed significant errors and 
instability in peak area calculations. 

3. The novel method could be applied to the activity measurement of LILW, and 
the reconstructed activity of Cs-137 was accurate, with an error of 2.4%, while 
the minimum measurement error of the Canberra system was 12.9%. 

4. The novel method could quantitatively analyze all of the nuclides in LILW 
without a nuclide library. For overlapping peaks composed of peak area ratios 
of 5:1:10, 10:1:20, and 20:1:40, the nuclide could be accurately identified; the 
area calculation error of the strong peaks was less than 1%, and the error of the 
weak peaks was less than 4%. 


This novel method can improve the accuracy and precision of LILW 
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measurements, provide key technical support for the reasonable disposal of LILW, 


and ensure the safety of humans and the environment. 


Acknowledgments 


The authors gratefully acknowledge the financial support provided by the 
National Natural Science Foundation of China (project number 12205190 and 
11805121) and the Science and Technology Commission of Shanghai Municipality 
(project number 21ZR1435400). 


References 


1. INTERNATIONAL ATOMIC ENERGY AGENCY, Classification of Radioactive 
Waste, IAEA Safety Standards Series No. GSG-1, IAEA, Vienna (2009) 


2. INTERNATIONAL ATOMIC ENERGY AGENCY, Interim Storage of 
Radioactive Waste Packages, Technical Reports Series No. 390, IAEA, Vienna 
(1998) 


3. INTERNATIONAL ATOMIC ENERGY AGENCY, Strategy and Methodology 
for Radioactive Waste Characterization, IAEA-TECDOC-1537, IAEA, Vienna 
(2007) 


4. P. Ormai, A. Fritz, J. Solymosi, et al., Inventory determination of low and 
intermediate level radioactive waste of Paks Nuclear Power Plant origin. J. 
Radioanal. Nucl. Chem. Artic. 211, 443-451 (1996). doi:10.1007/BF02039710 


5. M. Kashiwagi, H. Masui, Y. Denda, et al., ISO Standardization of the Scaling 
Factor Method for Low- and Intermediate Level Radioactive Wastes Generated at 
Nuclear Power Plants. In: //th International Conference on Environmental 
Remediation and Radioactive Waste Management, Parts A and B. ASMEDC; 
2007:625-629. doi:10.1115/ICEM2007-7015 


6. Standard Guide for Making Quality Nondestructive Assay Measurements, 
Accessed October 12, 2020. https://compass.astm.org/download/C1592C1592M- 
WITHDRAWN.4230.pdf 


7. American Society for Testing and Materials, Standard Test Method for 
Nondestructive Assay of Special Nuclear Material in Low Density Scrap and 
Waste by Segmented Passive Gamma-Ray Scanning. doi: 

25/31 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


I7: 


18. 


10.1520/C1133_C1133M-10 


Y.F. Bai, E. Mauerhofer, D.Z. Wang, et al., An improved method for the non- 
destructive characterization of radioactive waste by gamma scanning. Appl. 
Radiat. Isot. 67, 1897-1903 (2009). doi:10.1016/j.apradiso.2009.05.017 


Z.L. Cheng, Z. Lu, Z. Huang, et al., An improved segmented gamma scanning for 
radioactive waste drums. Math. Biosci. Eng. 16, 4491-4505 (2010). 
doi:10.3934/mbe.2019224 


Y.F. Bai, D.Z. Wang, M. Eric, et al., MC simulation of thermal neutron flux of 
large samples irradiated by 14 MeV neutrons. Nucl. Sci. Tech. 21(1):11-15(2010). 


C. Liu, W.G. Gu, N. Qian, et al., Study of image reconstruction using dynamic 
grids in tomographic gamma scanning. Nucl. Sci. Tech. 23,277—283 (2012) 


N. Qian, T. Krings, E. Mauerhofer, et al., Analytical calculation of the collimated 
detector response for the characterization of nuclear waste drums by segmented 
gamma scanning. J. Radioanal. Nucl. Chem. 292, 1325-1328 (2012). 
doi:10.1007/s10967-011-1601-1 


W.G. Gu, C. Liu, N. Qian, et al., Study on detection simplification of 
tomographic gamma scanning using dynamic grids applied in the emission 
reconstruction. Ann. Nucl. Energy. 58, 113-123 (2013). 
doi:10.1016/j.anucene.2013.02.031 


Q. Nan, D.Z. Wang, C. Wang, et al., Collimated LaBr3 detector response function 
in radioactivity analysis of nuclear waste drums. Nucl. Sci. Tech. 
24 ,060203(2013) 


W. Gu, K. Rao, D. Wang, et al., Semi-Tomographic Gamma Scanning Technique 
for Non-Destructive Assay of Radioactive Waste Drums. JEEE Trans. Nucl. Sci. 
63, 2793-2800 (2016). doi:10.1109/TNS.2016.2614964 


W.G. Gu, D.Z. Wang, X.H. Tang, et al., An Improved Gamma Scanning Assay 
Method for the 400-L Compacted Radioactive Waste Drum Based on the 
Segmented Equivalent Ring Source. JEEE Trans. Nucl. Sci. 66, 1889-1896 (2019). 
doi:10.1109/TNS.2019.2920180 


H. Yang, B. Dong, W.G. Gu et al., Transmission reconstruction algorithm by 
combining maximum-likelihood expectation maximization and a convolutional 
neural network for radioactive drum characterization. Appl. Radiat. Isot. 184, 
110172 (2022). doi:10.1016/j.apradiso.2022.110172 


H. Yang, H. Zhou, B. Dong, et al., A Novel Transmission Reconstruction 
Algorithm for Radioactive Drum Characterization. In: Volume 9: 


26/31 


19. 


20. 


21; 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


Decontamination and Decommissioning, Radiation Protection, and Waste 
Management. American Society of Mechanical Engineers; 2022:V009T09A005. 
doi:10.1115/ICONE29-90126 


Office of the Federal Register, National Archives and Records Administration. 10 
CFR 61 - LICENSING REQUIREMENTS FOR LAND DISPOSAL OF 
RADIOACTIVE WASTE. Code of Federal Regulations, title 10 :153-179(2016), 
https://www.govinfo.gov/app/details/CFR-2006-title 1 0-vol2/CFR-2006-title10- 
vol2-part61 Accessed 26 June 2022 


INTERNATIONAL ATOMIC ENERGY AGENCY, Determination and Use of 
Scaling Factors for Waste Characterization in Nuclear Power Plants, IAEA 
Nuclear Energy Series No. NW-T-1.18, IAEA, Vienna (2009) 


M.I. Abbas, M.S. Badawi, I.N. Ruskov, et al., Calibration of a single hexagonal 
Nal(TI) detector using a new numerical method based on the efficiency transfer 
method. Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. 
Assoc. Equip. 771, 110-114 (2015). doi:10.1016/j.nima.2014.10.056 


H. Yang, W. Gu, X. Zhang, et al., Research on the y spectrum-unfolding method 
of low- and intermediate-level radioactive waste based on LaBr3(Ce) detector. 
Radiat. Phys. Chem. 110841 (2023). doi:10.1016/j.radphyschem.2023.110841 


H. Yang, W. Gu, X. Zhang, et al., Research on the CdZnTe y spectrum analysis 
method based on an intelligent dynamic library. J. Radioanal. Nucl. Chem. 
(2023).Doi:10.1007/s10967-023-08858-9 


M.W. Rawool-Sullivan, R.J. Estep, D. Miko, Estimation of obliquely scattered 
gamma-ray response functions in the gross-count tomographic gamma scanner 
(GC-TGS) method. In: 1997 TEEE Nuclear Science Symposium Conference 
Record. IEEE; 1997:774-778. doi:10.1109/NSSMIC.1997.672697 


R. Venkataraman, M. Villani, S. Croft, et al., An integrated Tomographic Gamma 
Scanning system for non-destructive assay of radioactive waste. Nucl. Instrum. 
Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip. 579, 375- 
379 (2007). doi:10.1016/j.nima.2007.04.125 


J.T. Routti, S.G. Prussin, Photopeak method for the computer analysis of gamma- 
ray spectra from semiconductor detectors. Nucl. Instrum. Methods. 72, 125-142 
(1969). doi:10.1016/0029-554X(69)90148-7 


H. Cheng, Bao Hua Sun, L.H. Zhu, et al., Intrinsic background radiation of 
LaBr3(Ce) detector via coincidence measurements and simulations. Nucl. Sci. 
Tech. 31, 99 (2020). doi:10.1007/s41365-020-008 12-8 


M.A. Mariscotti, A method for automatic identification of peaks in the presence 


27/31 


29. 


30. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


of background and its application to spectrum analysis. Nucl. Instrum. Methods. 
50, 309-320 (1967). doi:10.1016/0029-554X(67)90058-4 


INTERNATIONAL ATOMIC ENERGY AGENCY, Intercomparison of Gamma 
Ray Analysis Software Packages, IAEA-TECDOC-1011, IAEA, Vienna (1998) 


Canberra. Genie 2000 Operations Manual. 
https://www3.nd.edu/~wzech/Genie%202000%200Operations%20Manual.pdf. 
Accessed 26 June 2022 


I.A. Slavic, S.P. Bingulac, A simple method for full automatic gamma-ray spectra 
analysis. Nucl. Instrum. Methods. 84, 261-268 (1970). doi:10.1016/0029- 
554X(70)90270-3 


A. Robertson, W.V. Prestwich, T.J. Kennett, An automatic peak-extraction 
technique. Nucl. Instrum. Methods. 100, 317-324 (1972). doi:10.1016/0029- 
554X(72)90701-X 


I.A. Slavić, Automatic analysis of gamma-ray spectra. Nucl. Instrum. Methods. 
112, 253-260 (1973). doi: 10.1016/0029-554X(73)90804-5 


H.P. Blok, J.C. de Lange, J.W. Schotman, A new peak search method for an 
automatic spectrum analysis program. Nucl. Instrum. Methods. 128, 545-556 
(1975). doi:10.1016/0029-554X(75)90523-6 


C.J. Sullivan, S.E. Garner, K.B. Butterfield, Wavelet analysis of gamma-ray 
spectra. In: JEEE Symposium Conference Record Nuclear Science 2004. Vol 1. 
IEEE; 2004:281-286. doi:10.1109/NSSMIC.2004.1462198 


J.W. Wang, W.G. Gu, H. Yang, D.Z. Wang. Analytical method for y energy 
spectrum of radioactive waste drum based on deep neural network. Nuclear 
techniques, vol. 45, no. 4, pp. 53-59(2022). 


ORTEC. GammaVision Users Manual. https://manualzz.com/doc/o/nbci8/ortec 
gammavision-v7-user-manual-783620h-mda--preset. Accessed 18 June 2022 


Clayton, ŒE. PIXAN: the Lucas Heights PIXE analysis computer 
package. AAEC/M--113(1986). ISBN 0 642 59841 X 


Miroslav Morhaé, An algorithm for determination of peak regions and baseline 
elimination in spectroscopic data. Nucl. Instrum. Methods Phys. Res. Sect. Accel. 
Spectrometers Detect. Assoc. Equip. 600, 478-487 (2009). 
doi:10.1016/j.nima.2008.11.132 


M.C. Ali Santoro, M.J. Anagnostakis, T. Boshkova, et al., Determining the 
probability of locating peaks using computerized peak-location methods in 


gamma-ray spectra as a function of the relative peak-area uncertainty. Appl. 
28/31 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


48. 


49. 


50. 


Sh 


Radiat. Isot. 155, 108920 (2020). doi:10.1016/j.apradiso.2019.108920 


Knoll, G.F. Radiation Detection and Measurement, 4th Edition, (Wiley, 
Hoboken,2010) 


International Electrotechnical (Commission. IEC 61453:2007 Nuclear 
instrumentation - Scintillation gamma ray detector systems for the assay of 
radionuclides - Calibration and routine tests. 
https://webstore.iec.ch/publication/5471. Accessed December 7, 2022 


M. Morháč, J.Kliman, V. Matoušek, et al., Efficient one- and two-dimensional 
gold deconvolution and its application to y-ray spectra decomposition. Nucl. 
Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip. 401, 
385-408 (1997). doi:10.1016/S0168-9002(97)01058-9 


C.M. Salgado, L.E.B. Brandão, R. Schirru, et al., Validation of a Nal(Tl) 
detector’s model developed with MCNP-X code. Prog. Nucl. Energy. 59, 19-25 
(2012). doi:10.1016/j.pnucene.2012.03.006 


M. Morháč, V. Matoušek, J. Kliman, Optimized multidimensional nonoscillating 
deconvolution. J. Comput. Appl. Math. 140, 639-658 (2002). doi:10.1016/S0377- 
0427(01)00521-0 


M. Morháč, Deconvolution methods and their applications in the analysis of -ray 
spectra. Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. 
Assoc. Equip. 559, 119-123 (2006). doi:10.1016/j.nima.2005.11.129 


Gold, R.AN ITERATIVE UNFOLDING METHOD FOR RESPONSE 
MATRICES. United States(1964).do1:10.2172/4634295. 


M.Morhaé, V.MatouSek, High-resolution boosted deconvolution of spectroscopic 
data. J: Comput. Appl. Math. 235, 1629-1640 (2011). 
doi:10.1016/j.cam.2010.09.005 


S. Steenstrup, A simple procedure for fitting a backgound to a certain class of 
measured spectra. J. Appl. Crystallogr. 14, 226-229 (1981). 
doi:10.1107/S002188988 1009242 


S.A. Gerasimov, Simple method for background subtraction in gamma-ray 
spectra. Int. J. Rad. Appl. Instrum. [A]. 43, 1529-1531 (1992). doi:10.1016/0883- 
2889(92)90186-I 


C.V. Hampton, B. Lian, W.C. McHarris. Fast-Fourier-transform spectral 
enhancement techniques for y-ray spectroscopy. Nuclear Instruments and 
Methods in Physics Research Section A, Accelerators, Spectrometers, Detectors 
and Associated Equipment, 353(1-3), 280-284(1994). ISSN 0168-9002 


29/31 


52. 


53. 


54. 


55. 


56. 


57. 


58. 


59. 


60. 


61. 


62. 


63. 


C.G. Ryan, E. Clayton, W.L. Griffin, et al., SNIP, a statistics-sensitive 
background treatment for the quantitative analysis of PIXE spectra in geoscience 
applications. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 
34, 396-402 (1988). doi:10.1016/0168-583X(88)90063-8 


C.J. Sullivan, M.E. Martinez, S.E. Garner. Wavelet Analysis of Sodium Iodide 
Spectra. IEEE Nucl. Sci. Symp. Conf. Rec. Wyndham El Conquistador Resort, 
Puerto Rico: IEEE; p. 302—6(2005). doi:10.1109/NSSMIC.2005.1596258. 


S. Lovejoy, D. Schertzer, Haar wavelets, fluctuations and structure functions: 
convenient choices for geophysics. Nonlinear Process. Geophys. 19, 513-527 
(2012). doi:10.5194/npg-19-513-2012 


J.M. Gregoire, D. Dale, R.B. van Dover, A wavelet transform algorithm for peak 
detection and application to powder x-ray diffraction data. Rev. Sci. Instrum. 82, 
015105 (2011). doi:10.1063/1.3505103 


J. Uher, G. Roach, J. Tickner, Peak fitting and identification software library for 
high resolution gamma-ray spectra. Nucl. Instrum. Methods Phys. Res. Sect. 
Accel. Spectrometers Detect. Assoc. Equip. 619, 457-459 (2010). 
doi:10.1016/j.nima.2009.12.086 


A. Likar, A peak-search method based on spectrum convolution. J. Phys. Appl. 
Phys. 36, 1903-1909 (2003). doi:10.1088/0022-3727/36/15/323 


M. Morhaé, J. Kliman, V. MatouSek, et al., Background elimination methods for 
multidimensional coincidence y-ray spectra. Nucl. Instrum. Methods Phys. Res. 
Sect. Accel. Spectrometers Detect. Assoc. Equip. 401, 113-132 (1997). 
doi:10.1016/S0168-9002(97)01023-1 


J.L. Coolidge, The Story of the Binomial Theorem. Am. Math. Mon. 56, 147-157 
(1949). doi:10.1080/00029890.1949.11999350 


M.S. Badawi, A.M. El-Khatib, M.E. Krar, New numerical simulation approach to 
calibrate the Nal(T1) detectors array using non-axial extended spherical sources. J. 
Instrum. 8, P11005-P11005 (2013). doi:10.1088/1748-0221/8/11/P11005 


M.S. Badawi, M. Abd-Elzaher, A.A. Thabet, et al., An empirical formula to 
calculate the full energy peak efficiency of scintillation detectors. Appl. Radiat. 
Isot. 74, 46-49 (2013). doi:10.1016/.apradiso.2012.12.011 


Canberra.NDA 2000. https://www.mirion.com/products/nda-2000-non- 
destructive-assay-software. Accessed 15 Nov.2022 


American Society for Testing and Materials. ASTM E181-17 Standard Test 
Methods for Detector Calibration and Analysis of Radionuclides. doi: 


30/31 


64. 


65. 


66. 


67. 


68. 


69. 


70. 


71. 


72. 


73. 


10.1520/E0181-17 


ORTEC. Comparison of Gamma Ray Nondestructive Assay Measurement 
Techniques. https://www.ortec-online.com/-/media/ametekortec/technical- 
papers/waste-assay/comparison-gamma-ray-nondestructive-assay-measurement- 
techniques.pdf?la=en&revision=4 | 5c3b6e-e8 13-4930-abb1- 
c79c5870c644&hash=2D 10984B0F 16B476B43B17F6A332258B. Accessed 19 
Nov.2022 


S. Agostinelli, J. Allison, K. Amako, et al., Geant4—a simulation toolkit. Nucl. 
Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip. 506, 
250-303 (2003). doi:10.1016/S0168-9002(03)01368-8 


W. Tang, J.G. Liang, Y. Ge, et al., A method for neutron-induced gamma spectra 
decomposition analysis based on Geant4 simulation. Nucl. Sci. Tech. 33, 154 
(2022). doi:10.1007/s41365-022-01144-5 


A.M. El-Khatib, M.S. Badawi, A.A. Thabet, et al., Well-type Nal(TI) detector 
efficiency using analytical technique and ANGLE 4 software based on 
radioactive point sources located out the well cavity. Chin. J. Phys. 54, 338-346 
(2016). doi:10.1016/j.cjph.2016.03.019 


M.T. Haj-Heidari, M.J. Safari, H. Afarideh, et al., Method for developing HPGe 
detector model in Monte Carlo simulation codes. Radiat. Meas. 88, 1-6 (2016). 
doi:10.1016/j.radmeas.2016.02.035 


Y.M. Gomez, H.A.C. Aguila, C.M.A. Hernandez, et al., Monte Carlo simulation 
of the efficiency response of a well-type HPGe detector at 46.54 keV (2015). 


http://scielo.sld.cu/pdf/nuc/n57/nuc015715.pdf .Accessed 28 Jan.2023 


M.S. Badawi, S.I. Jovanovic, A.A. Thabet, et al., Calibration of 42 Nal(T1) 
detectors with coincidence summing correction using new numerical procedure 
and ANGLE4 software. AIP Adv. 7, 035005 (2017). doi:10.1063/1.4978214 


M. JeSkovsky, A. Javorník, R. Breier, et al., Experimental and Monte Carlo 
determination of HPGe detector efficiency. J. Radioanal. Nucl. Chem. 322, 1863- 
1869 (2019). doi:10.1007/s10967-019-06856-4 


American Society for Testing and Materials. American National Standard 
Calibration and Use of Germanium Spectrometers for the Measurement of 
Gamma-Ray Emission Rates of Radionuclides. ANSI Std N42.14-1991, vol., no., 
pp. 1-76, 30(1991). doi: 10.1109/IEEESTD.1991.101072. 


IEEE Standard Test Procedures for Germanium Gamma-Ray Detectors, JEEE Std 
325-1996. 0_1- Published online 1997. doi:10.1109/IEEESTD.1997.82400 


31/31 


